%%%%%%%1-HIDROXYPYRENE
X(1,:,:)=(etohp);
X(2,:,:)=(pyr101ppm);
X(3,:,:)=(pyr1001ppm);
X(4,:,:)=(pyr10001ppm);
X(5,:,:)=(pyr100001ppm);
X(6,:,:)=(pyr1000001ppm);

Xin.X = X;
Xin.Em = em;
Xin.Ex = ex;

EEMcor=smootheem(Xin,[20 10],[5 15],[20 15],[10 15],[],[],0.5);

EEMP1(1,:,:)=EEMcor.X(2,:,:)-EEMcor.X(1,:,:);
EEMP1(2,:,:)=EEMcor.X(3,:,:)-EEMcor.X(1,:,:);
EEMP1(3,:,:)=EEMcor.X(4,:,:)-EEMcor.X(1,:,:);
EEMP1(4,:,:)=EEMcor.X(5,:,:)-EEMcor.X(1,:,:);
EEMP1(5,:,:)=EEMcor.X(6,:,:)-EEMcor.X(1,:,:);

[ssX,Corco,It] = pftest(5,EEMP1,5,[],[2 2 2])

%%%%%%%1 Y 2-HIDROXYNAPHTHALENE
X(1,:,:)=(etohp);
X(2,:,:)=(naf101ppm);
X(3,:,:)=(naf10075ppm);
X(4,:,:)=(naf1001ppm);
X(5,:,:)=(naf10001ppm);
X(6,:,:)=(naf100001ppm);
X(7,:,:)=(naf201ppm);
X(8,:,:)=(naf2001ppm);
X(9,:,:)=(naf20001ppm);
X(10,:,:)=(naf200001ppm);
X(11,:,:)=(naf2000001ppm);

Xin.X = X;
Xin.Em = em;
Xin.Ex = ex;

EEMcor=smootheem(Xin,[20 10],[5 15],[20 15],[10 15],[],[],0.5);

EEMN2(1,:,:)=EEMcor.X(2,:,:)-EEMcor.X(1,:,:);
EEMN2(2,:,:)=EEMcor.X(3,:,:)-EEMcor.X(1,:,:);
EEMN2(3,:,:)=EEMcor.X(4,:,:)-EEMcor.X(1,:,:);
EEMN2(4,:,:)=EEMcor.X(5,:,:)-EEMcor.X(1,:,:);
EEMN2(5,:,:)=EEMcor.X(6,:,:)-EEMcor.X(1,:,:);
EEMN2(6,:,:)=EEMcor.X(7,:,:)-EEMcor.X(1,:,:);
EEMN2(7,:,:)=EEMcor.X(8,:,:)-EEMcor.X(1,:,:);
EEMN2(8,:,:)=EEMcor.X(9,:,:)-EEMcor.X(1,:,:);
EEMN2(9,:,:)=EEMcor.X(10,:,:)-EEMcor.X(1,:,:);
EEMN2(10,:,:)=EEMcor.X(11,:,:)-EEMcor.X(1,:,:);

[ssX,Corco,It] = pftest(10,EEMN2,5,[],[2 2 2])

EEMN(1,:,:)=EEMcor.X(2,:,:)-EEMcor.X(1,:,:);
EEMN(2,:,:)=EEMcor.X(3,:,:)-EEMcor.X(1,:,:);
EEMN(3,:,:)=EEMcor.X(4,:,:)-EEMcor.X(1,:,:);
EEMN(4,:,:)=EEMcor.X(5,:,:)-EEMcor.X(1,:,:);
EEMN(5,:,:)=EEMcor.X(6,:,:)-EEMcor.X(1,:,:);
EEMN(6,:,:)=EEMcor.X(7,:,:)-EEMcor.X(1,:,:);
EEMN(7,:,:)=EEMcor.X(8,:,:)-EEMcor.X(1,:,:);
EEMN(8,:,:)=EEMcor.X(9,:,:)-EEMcor.X(1,:,:);
EEMN(9,:,:)=EEMcor.X(10,:,:)-EEMcor.X(1,:,:);
EEMN(10,:,:)=EEMcor.X(11:,:)-EEMcor.X(1,:,:);

[ssX,Corco,It] = pftest(9,EEMN,5,[],[2 2 2])

%%%%%%%9-HIDROXYPHENATHRENE
X(1,:,:)=(etohp);
X(2,:,:)=(fen91ppm);
X(3,:,:)=(fen901ppm);
X(4,:,:)=(fen9001ppm);
X(5,:,:)=(fen90001ppm);

Xin.X = X;
Xin.Em = em;
Xin.Ex = ex;

EEMcor=smootheem(Xin,[20 10],[5 15],[20 15],[10 15],[],[],0.5);

EEMPH(1,:,:)=EEMcor.X(2,:,:)-EEMcor.X(1,:,:);
EEMPH(2,:,:)=EEMcor.X(3,:,:)-EEMcor.X(1,:,:);
EEMPH(3,:,:)=EEMcor.X(4,:,:)-EEMcor.X(1,:,:);
EEMPH(4,:,:)=EEMcor.X(5,:,:)-EEMcor.X(1,:,:);

[ssX,Corco,It] = pftest(4,EEMPH,5,[],[2 2 2])
